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ABSTRACT 

Core collapse of massive stars with a relativistic jet expulsion along the rotation 
axis is a widely discussed scenario for gamma-ray burst (GRB) production. How- 
ever the nature of the stellar progenitor remains unclear We study the evolution of 
stars that may be the progenitors of long-soft GRBs - rotating naked helium stars 
presumed to have lost their envelopes to winds or companions. Our aim is to in- 
vestigate the formation and development of single and binary systems and from this 
population evaluate the rates of interesting individual species. Using a rapid binary 
evolution algorithm, that enables us to model the most complex binary systems and 
to explore the effect of metallicity on GRB production, we draw the following con- 
clusions. First we find that, if we include an approximate treatment of angular mo- 
mentum transport by mass loss, the resulting spin rates for single stars become too 
low to form a centrifugally supported disc that can drive a GRB engine - although 
they do have sufficiently massive cores to form black holes. Second massive stars 
in binaries result in enough angular momentum - due to spin-orbit tidal interactions 
- to form a centrifugally supported disc and are thus capable of supplying a suffi- 
cient number of progenitors. This holds true even if only a small fraction of bursts 
are visible to a given observer and the GRB rate is several hundred times larger than 
the observed rate. Third low-metallicity stars aid the formation of a rapidly rotating, 
massive helium cores at collapse and so their evolution is likely to be affected by 
the local properties of the ISM. This effect could increase the GRB formation rate 
by a factor of 5-7 al Z = Zq/2{){). Finally we quantify the effects of mass loss, 
common-envelope evolution and black-hole formation and show that more stringent 
constraints to many of these evolution parameters are needed in order to draw quan- 
titative conclusions from population synthesis work. 
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1 INTRODUCTION 



Gamma-ray bursts (GRBs) are intense flashes of gamma-rays which, for a few seconds, light up 
an otherwise dark gamma-ray sky. They are detected about once a day and outshine every other 
gamma-ray source in the sky. Major advances have been made in the last few years including 
the discovery of slowly fading X-ray (Costa et al. 1997), optical (van Paradijs et al. 1997) and 
radio (Frail et al. 19997) afterglows of GRBs, the identification of host galaxies at cosmological 
distances and the evidence that many of them are associated with star forming regions and possibly 
supemovae (see Meszaros 2001 for a recent review). Given the twin requirements of enormous 
energy, about 10^'^ erg, and association with star forming regions, the currently favoured models 
all involve massive, collapsing stars and their end products, especially black holes. 

Hyper-accreting black holes (0.01 to 10 Mq s~^) can arise from coUapsars of various types 
(Woosley 1993; Fryer et al. 2001); white dwarf-black hole mergers (Frail et al. 1999; Belczynski 
et al. 2002a); helium core-black hole mergers (Fryer & Woosley 1998; Belczynski et al. 2002a); 
and merging neutron stars and black holes (Kluzniak & Lee 1998; Ruffert & Janka 1999; Lee 
& Ramirez-Ruiz 2002). Those models that involve merging neutron stars and black holes mostly 
transpire outside star forming regions and are not currently thought to be appropriate for long 
GRBs (Meszaros 2001). Of the other accretion-based models, black holes merging with both he- 
lium cores and white dwarfs have the advantage that they could occur in star forming regions 
(Belczynski et al. 2002a), produce adequate energy, and have sufficient angular momentum to 
form a disc. However, if the burst is to be beamed to less than 1% of the sky, the white dwarf 
model may provide an inadequate number of events (Fryer et al. 1999). The helium core model 
requires that a hydrogen envelope initially be present in order that the black hole (or neutron star) 
experience a common envelope evolution and merge with the core. Yet the same envelope must be 
absent when the burst occurs or prohibitive baryon loading will occur (MacFadyen et al. 2001). 
It may be that the merger peels off the envelope as the compact remnant goes in, but calculations 
to support this hypothesis are needed. It should be noted that the merger of a black hole with a 
massive helium core produces essentially the same conditions as in ordinary coUapsars, except for 
a larger angular momentum. The beaming and energetics should also be comparable. 

A coUapsar forms when the evolved core of a massive star collapses to a black hole, either by 
fallback or because the iron core fails to produce an outgoing shock (Woosley 1993; Paczyhski 
1998; MacFadyen & Woosley 1999). Prompt and delayed black hole formation may occur in stars 
with a range of radii depending on the evolutionary state of the massive progenitor, its metallicity 
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and its multiplicity. The two necessary ingredients for the coUapsar are a failed or weak initial 
supernova (SN) explosion that produces a black hole and sufficient angular momentum to form a 
disc. Energy dissipated in the disc or the rotation of the black hole itself is assumed to power a 
jet of high Lorentz factor F 100 (Soderberg & Ramirez-Ruiz 2002) which escapes the stellar 
progenitor along the polar axis. 

The traversal time for the relativistic jet through the hydrogen envelope of a typical massive 
star ranges from hundreds to thousands of seconds. At the time of the burst naked helium stars, 
which have radii of only a few light seconds, are required if the lifetime of the engine is not 
to be short when compared with the time it takes the jet to tunnel through the star (MacFadyen 
& Woosley 1999; Aloy et al. 2000; Wheeler et al. 2000; Meszaros & Rees 2001; Meszaros & 
Waxman 2001; Ramirez-Ruiz, Celotti & Rees 2002; Matzner 2003). This is aided if the stellar 
progenitor undergoes a Wolf-Rayet (WR) phase characterised by a strong stellar wind that causes 
the star to lose enough of its outer layers for the surface hydrogen abundance to become minimal 
(MacFadyen & Woosley 1999; Ramirez-Ruiz et al. 2001). The radius of a WR star is then suffi- 
ciently small for the jet to emerge before the engine ceases to operate. Break-out is further helped 
if the progenitor originates in a low-metallicity environment, when stars are smaller in radius, lose 
less mass and consequently have more massive cores. Both effects are thought to favour GRB 
formation (MacFadyen & Woosley 1999; Fryer et al. 1999). 

Here we present population synthesis calculations which allow us to explore what effect initial 
mass, metallicity and membership in a binary system have on GRB production. In addition to all 
aspects of single-star evolution, our rapid evolution algorithm includes features such as mass trans- 
fer, mass accretion, common-envelope evolution, collisions, SN kicks and angular-momentum- 
loss mechanisms which enables us to model even the most complex binary systems. The binary- 
evolution algorithm is described in detail in Section 2. Section 3 contains the results of binary and 
single populations of core collapse SNe in solar metallicity environments. We also quantify the 
effects that mass loss, common-envelope evolution, and black-hole formation have on GRB pro- 
duction. The various kinds of companion remnants that are expected when the primary massive 
star collapses to a black hole are also discussed in Section 3. We present the results of population 
synthesis to examine the effects of metallicity in Section 4 and compare our model results with 
observations in Section 5. Our conclusions are then given in section 6. 
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2 STELLAR EVOLUTION MODELS 



2.1 A field guide to stellar evolution 

A fundamental tool in stellar evolution is the Hertzsprung-Russell diagram (HRD) which provides 
a correspondence between the observable stellar properties - the effective surface temperature 
Teffand luminosity L. Fig. 1 shows the evolution of a collection of stars in the HRD from the zero- 
age main sequence (ZAMS), where a star adjusts itself to nuclear and hydrostatic equilibrium, 
until the end of nuclear burning. As stars take a relatively short time to reach the ZAMS after 
molecular cloud collapse, all ages are measured from this point. The span of a star's life, its path 
on the HRD and its ultimate fate depend critically on its mass. 

The immediate post-MS evolution towards the right in the HRD occurs at nearly constant 
luminosity and is very rapid. For this reason very few stars are seen in this phase, and this region 
of the HRD is called the Hertzsprung gap (HG) or the sub-giant branch. During this HG phase the 
radius of the star increases greatly causing a decrease in Tgfj. As the convective envelope grows in 
extent the star reaches the giant branch (GB). Eventually a point is reached on the GB where the 
core temperature is high enough for stars to ignite their central helium supply. When core helium 
burning (CHeB) begins the star descends along the GB until contraction moves the star away from 
the fully convective region of the HRD and back towards the MS in what is called a blue loop. 
During CHeB, carbon and oxygen are produced in the core. Eventually core helium is exhausted 
and the star moves back to the right in the HRD. Stars of very high mass (M^ams > SSM©) reach 
high enough central temperatures on the HG for helium to ignite before reaching the GB. 

Evolution after the exhaustion of core helium is very similar to that occurring after core- 
hydrogen exhaustion at the end of the MS - the convective envelope deepens again to begin what 
is called the asymptotic giant branch (AGB). On the AGB the star is now a dense core composed 
of carbon and oxygen surrounded by a helium burning shell which in turn adds carbon (and some 
oxygen) to the degenerate core. Initially the H-buming shell is weak so that the luminosity is 
supplied by the He-burning shell. During this phase, the stellar radius grows significantly. 

The surface gravity of the star is lowered so the surface material is less tightly bound. Mass 
loss from the stellar surface becomes significant with the rate of mass loss increasing with time. 
Unfortunately our understanding of the mechanisms that cause this mass loss is poor with possible 
suggestions linking it to the helium shell flashes or to periodic envelope pulsations. Whatever the 
cause, the influence on the evolution of AGB stars is significant. Mass loss eventually removes all 
of the star's envelope so that the H-buming shell shines through. 
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If the mass of the star is large enough, Mzams ^ ^Mq (although the exact value depends on 
its metallicity and mass-loss history), the carbon-oxygen core is not very degenerate and carbon 
ignites as it contracts, followed by a succession of nuclear reaction sequences which very quickly 
produce an inner iron core. Any further reactions are endothermic and cannot contribute to the 
luminosity of the star. Photo-disintegration of iron, combined with electron capture by protons 
and heavy nuclei, then removes most of the electron degeneracy pressure supporting the core and 
it begins to collapse rapidly. When the density becomes large enough the inner core rebounds 
sending a shock wave outwards through the outer layers of the star that have remained suspended 
above the collapsing core. As a result the envelope of the star is ejected in a SN explosion so 
that the AGB is truncated soon after the start of carbon burning. The remnant in the inner core 
stabilises to form a neutron star (NS) supported by neutron degeneracy pressure unless its mass is 
large enough that complete collapse to a black hole occurs. 

Stars with Afzams ^ 25M0 are severely affected by mass loss during their entire evolution and 
may lose their hydrogen envelopes, exposing nuclear processed material. If this occurs a naked 
helium star is produced and such stars (or stars about to become naked helium stars) may be 
WR stars. WR stars are massive objects which are found near the MS, are losing mass at very high 
rates and show weak or no hydrogen lines in their spectra. Naked helium stars can also be produced 
from less massive stars in binaries as a consequence of mass transfer (see Fig. lb). Variation in 
composition affects the stellar evolution timescales as well as the appearance in the HRD and even 
the ultimate fate of the star. A more detailed discussion of the various phases of evolution can be 
found in Hurley at el. (2000). 

2.2 Binary evolution 

The evolution of binary stars does not differ from that of single stars unless they get in each other's 
way. If the binary orbit is wide enough the individual stars are not affected by the presence of a 
companion so that standard stellar evolution theory is all that is required to describe their devel- 
opment. However if the stars become close they can interact, with consequences for the evolution. 
Models for binary evolution have been implemented in the past (e.g. Whyte & Eggleton 1985; Pols 
& Marinus 1994; Portegies Zwart & Verbunt 1996). Here we carry out computational simulations 
with the rapid-binary evolution algorithm first developed by Tout et al. (1997) and recently updated 
by Hurley et al. (2002). Amongst other improvements Hurley et al. (2002) incorporate the detailed 
single-star evolution formulae of Hurley, Pols & Tout (2000) which allow for a wider range of 
stellar types than the description of stellar evolution used by Tout et al. (1997). This requires an 
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update of the treatment of processes such as Roche lobe overflow, common-envelope evolution 
and coalescence by collision (the reader is referred to Hurley et al. 2002 for further details). 

As a result of Roche-lobe overflow it is possible for the binary components to come into 
contact and coalesce or for the binary to reach a common-envelope (CE) state. The most frequent 
case of common-envelope evolution involves a giant transferring mass to a main-sequence star on 
a dynamical timescale. Although the process is difficult to model, and therefore uncertain, it is 
envisaged that the secondary is not able to accept the overflowing material owing to its relatively 
long thermal timescale. The giant envelope overfills the Roche-lobes of both stars so that the 
giant core and the MS star are contained within a common-envelope. Owing to its expansion the 
envelope rotates slower than the orbit of the core and the MS star so that friction causes them to 
spiral together and transfer energy to the envelope. We assume that the cores spiral-in, transferring 
orbital energy to the envelope with an efficiency oce, which is necessarily a free parameter due to 
uncertainty in its value. It is probably not a constant (Regos & Tout 1995) but generally ctcE ~ 1 
is used. 

Throughout this paper we refer to one star as the primary, mass Mi, and the other as the sec- 
ondary (or companion), mass M2. At any time the primary is the star filling, or closest to filling, 
its Roche lobe. Numerical values of mass, luminosity and radius are in solar units unless indicated 
otherwise. The algorithm provides the stellar luminosity, radius, core mass, core radius, and angu- 
lar momentum, for each of the component stars as they evolve. A prescription for mass loss from 
stellar winds is also included in the algorithm. The algorithm covers all the evolution phases from 
zero-age main sequence (ZAMS) up to and including the remnant stages and is valid for all masses 
in the range 0.1 to 100 Mq and metallicities from Z — 10~^ to 0.03. This rapid binary-evolution 
algorithm is a natural extension of the single-star evolution algorithm. Our aim is to evolve a pop- 
ulation of binaries according to chosen distributions of primary mass, secondary mass and orbital 
separation, in conjunction with a realistic birth rate function and, from this population, to calculate 
birth rates and expected numbers in the Galaxy for stars that explode as type Ib/c SNe, thought to 
be possible GRB progenitors. 



2.3 Method 

We first set up a grid of initial binary parameters (in the detailed single-star evolution the secondary 
is removed from the system) Mi, M2 and separation a within the limits: 

Mi,M2: O.l^SO.OMo (1) 
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a: 3.0^10^Rq (2) 
with the — 100 grid points of parameter x logarithmically spaced, 

5 In X = ^—r [In Xmax - In Xmin] ■ (3) 

For each set of initial parameters we evolve the binary system to an age of 100 Myr. If a binary 
system j evolves through a phase that is to be identified with a certain individual binary population 
i then the system makes a contribution 

Srj = 5$ (In Mij )(p{\n ) * (In % ) 5 In Mi 5 In Ms 5 In a (4) 

to the rate rj at which that particular population is born. This rate depends on the star formation 
rate S, the primary mass distribution $ (In Mi), the secondary mass distribution ip (InMs) and 
the separation distribution ^ (In a). A star formation rate of 5" = 7.608 galaxy"^ yr~^ is used in 
this work which corresponds to one star with M > 0.8 forming in the Galaxy every year 
(see Hurley et al. 2002). Such a rate is in rough agreement with the birth rate of white dwarfs 
in the Galaxy Ri 2 x 10"^^ pc'^ yr"^ (Phillips 1989), noting that only stars with M > 0.8 M© 
can possibly evolve to white dwarfs in the age of the Galaxy^, and assuming an effective Galactic 
volume of 5 X 10^^ pc'^ The primary mass distribution ^{m) is the IMF of Kroupa, Tout & 
Gilmore (1993; hereafter KTG), 

m ^ mo 

, , al■m~^^^ mo < m ^ 0.5 

eM = <! (5) 

a2'm 0.5 < m ^ 1.0 
a2m~^^'^ 1.0 < m < oo, 
where ^ (m) dm is the probability that a star has a mass, expressed in solar units, between m and 
m + dm. The distribution is normalized according to ^{m) dm — 1, so that, for mo = 0.1, 
oi = 0.29056 and = 0.15571. Then $ (In Mi) = Mi ^ (Mi). If the component masses are to be 
chosen independently, the secondary mass distribution is then (/9 (In M2) = M2 ^ (M2). However 
there is observational evidence (Eggleton, Fitchett & Tout 1989; Mazeh et al. 1992; Goldberg 

6 Mazeh 1994) to support correlated masses: ip (InMs) = ^ = q2, which corresponds to a 
uniform distribution of the mass-ratio q2, for < ^2 ^ 1 (the flat-q distribution). The separation 
distribution is taken to be ^ (In a) — '&,a constant, between the limits 3 and 10^ it!©. Normalization 
gives 1? = 0.12328. 



^ The same assumptions regarding the star formation rate have been made previously (e.g. Iben & Tutukov 1984; Han 1998; Hurley et al. 2002) 
facilitating comparison with these results. 
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2.4 Mass loss and rotation 

The particular mass-loss prescriptions through the various phases, which have been found to fit 
the observations well, are described in detail by Hurley et al. (2000; 2002). On the GB phase 
and beyond we apply mass loss to the envelope according to the formula of Kudritzki & Reimers 

(1978), while for the AGB we apply the formulation of Vassiliadis & Wood (1993). For massive 
stars we model mass loss over the entire HRD using the prescription given by Nieuwenhuijzen & 
de Jager(1990), 

M^j = (^I^^ 9.6 X 10-^^i?°-«iLi-2^M°-^6 Meyr-i (6) 

for L > 4OOOL0, modified by the factor Z^/^ (Kudritzki et al. 1989). Numerical values of mass, 
luminosity and radius are in solar units, unless otherwise specified. 

For small H-envelope mass, fi < 1.0, we include a Wolf-Rayet-like mass loss (Hamann & 
Koesterke 1998) which we have reduced to give 

MwR = lO-'^L'-' (1.0 - 1^) Moyr-i (7) 

where fx is given by 
Menv 



5.0,maxn.2, 



(8) 



L 

To. 

Here Lq = 7.0 x 10^ Lq and k — —0.5 (Hurley et al. 2000). The reduction in Mwr is necessary 
in order to produce sufficient black holes to match the number observed in binaries (see Hurley et 
al. 2000). 

As we plan to use the evolution routines to investigate GRB production, it is desirable to follow 

the stars' angular momentum. To do this we must start each star with a realistic spin on the ZAMS. 

A reasonable fit to the v^ot MS data of Lang (1992) is given by 
330M^-^ 

'-"(^^= 15.0 + M3.45 ^"^^" (9) 
so that 

Q = 45.35-^^ yr"^ . (10) 

^zams 

The angular momentum is then given by 

J,pi^^ m^kMR^n (11) 

where the constant k depends on the internal structure, e.g. k = 2/5 for a solid sphere and k = 2/3 
for a spherical shell. In actual fact we find the angular momentum by splitting the star into two 
parts, consisting of the core (with mass Mc) and the envelope (with mass Menv), so that 
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Jspin = {k2 (M - Mc) + hM,R,^) n (12) 

where k2 = 0.1, based on detailed giant models which reveal k = O.lMenv/M, and fca = 0.21 for 
ann = 3/2 polytrope such as a white dwarf, neutron star or dense convective core. This works well 
for post-MS stars which have developed a dense core whose rotation is likely to have decoupled 
from the envelope while also representing the near uniform rotation of homogeneous MS stars 
which have no core. When the star loses mass its stellar wind carries off angular momentum at a 
rate given by 

j = kMh (13) 
where h — E?Q,. Thus 

Jspin(t + At) = J,pin(i) - ^AMi?2Q (14) 

when the star loses an amount of mass AM, where we take A; = 2/3 because we assume that all 
the mass is lost uniformly at the surface of the star ie. from a spherical shell. 

We also include magnetic braking for stars that have appreciable convective envelopes 

j^b = 5.83 X 10-'"^ {mf MeR^yr-^, (15) 

with Q, in units of years. Following Rappaport et al. (1983) we don't allow magnetic braking for 
fully convective stars, M < O.SSM© but this restriction does not affect our results here. 



2.5 Core collapse progenitors 

Core-collapse SNe spectra may contain hydrogen absorption lines (type II) or show very weak or 
no hydrogen lines at all (type Vole) thought to indicate the lack of a hydrogen envelope around 
the imploding star. As in Hurley et al. (2002), we define a type Ib/c SN as one that produces 
a neutron star or black hole from a primary, naked helium star. All other SN progenitors that 
produce a neutron star or black hole are considered to be type II. The carbon-oxygen core mass 
Mqo of a SN Ib/c progenitor has a maximum value of maxfAfcin 0.773Mhc — 0.35 M©] (Hurley 
2000; hereafter HOO), where Mne is the initial mass of the helium star and Mch = 1.44 Mq is 
the Chandrasekhar mass. For masses above this, a SN explosion is assumed to take place. The 
expression for the maximum mass of the core of a SN II progenitor is similar with Mne replaced 
by the mass of the carbon-oxygen core at the base of the asymptotic giant branch. If 

Mco > 7.0 M© (16) 
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we assume a black hole forms (Woosley 1993). The mass of the neutron star (NS) or black hole 
(BH) formed is then given by (HOO) 

Mns/bh = 1.17 + O.OQMco Mq (17) 

corresponding to neutron star baryonic masses in the range 1.3 to 1.8 M© and a minimum 
black hole mass of 1.8 M0. 

A Wolf-Rayet stellar wind is used in the stellar models when the H-envelope mass 

Menv is small jji < 1.0. Stars with n > 1.0 are excluded from the Wolf-Rayet population. 

We apply the condition that a star which is to be considered as a GRB progenitor must: (i) 
go through a Wolf-Rayet phase, as defined above, before ending its life as a type Ib/c SN; (ii) 
have a sufficiently massive core to form a black hole; and (iii) have enough angular momentum at 
the time of collapse to allow the formation of a disc. To decide whether a centrifugally supported 
accretion disc can form around a central black hole we compare the angular momentum of the 
progenitor star to that which a test particle would require at the last stable orbit (LSO) around a 
black hole (see e.g. Heger & Woosley 2002). In the case of a non-rotating, spherically symmetric 
black hole the last stable orbit has a radius tlso = GGMbh/ c^. The condition that the matter does 
not spiral into the black hole is simply that the rotational energy of the matter should be greater or 
equal to half the gravitational potential energy: | J^// ^ |GMBH"^/rLso> where J is the angular 
momentum of the element mass m and / = nT'^tso is its moment of inertia about the hole in the 
LSO. In other words, j > GMbh?^lsO' where j is the specific angular momentum j — J/m. The 
inner boundary condition is therefore that the angular momentum with which the material arrives 
at the stable orbit should be greater than or equal to this critical value. 

3 PROGENITORS OF CORE COLLAPSE SNe AND GRBs 
3.1 Supernovae 

The effects of close binary evolution are observed in many systems, such as cataclysmic variables. 
X-ray binaries and Algols and by the presence of stars such as blue stragglers which cannot be ex- 
plained by single star evolution. The fraction of star systems containing at least two gravitationally 
bound stars, i.e. the binary fraction, is at least 50% (see Larson 2001 for a recent review). In the 

case of Wolf-Rayet stars it is perhaps as high as 80%. Some uncertainty still exists due to small 
number statistics but multiplicity seems a crucial element of massive star evolution. 

While many of the processes involved are not understood in detail we do have a qualitative 
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picture of how binaries evolve and can hope to construct a model that correctly follows them 
through the various phases of evolution. In order to conduct statistical studies of complete binary 
populations, i.e. population synthesis, such a model must be able to produce any type of binary 
that is observed in enough detail. 

For most of the individual binary populations, observational birth rates or numbers in the 
Galaxy are uncertain because of selection effects involved when undertaking surveys. However 
enough data exist overall to enable a meaningful comparison with the results. Cappellaro (2001) 
combined the results of five independent SN searches to obtain a sample of 137 SNe from which 
he derived the following birth rates for the Galaxy: 

2.4 ± 1 X 10-3 yr-i SNIb/c 
14.9 ± 6 X 10-3 yr-i SNH. 

The values reported above assume Hq = 65 km s^^ Mpc^\ and that the Galaxy has an average 
SN rate its morphological type (here assumed Sb-Sbc) and luminosity (2.3 x 10^°Lb see Cap- 
pellaro et al. 1999). Our predicted formation rates for SNe Ib/c and SNe II in single and binary 
stellar populations are shown in Table 1. Binary models have been evolved according to chosen 
distributions of primary (KTG distribution) and secondary (KTG or q-flat distribution) masses. 
Throughout this paper we refer to these distributions as KTG-KTG or KTG-q, respectively. The 
single star results are consistent with the observed rates, as are the binary type II SNe at Z = 0.02. 
A discrepancy lies in the birth rates of type Ib/c SNe because the theoretical rates in the KTG-q 
case are high but not for the KTG-KTG case which are in good agreement. This can be reconciled 
by assuming a binary fraction smaller than unity (Hurley et al. 2002). 

The SN II rate is a fairly robust number but the SN Ib/c rate is sensitive to many of the assump- 
tions underlying the model (Hurley et al. 2002). In binary models the majority of type Ib/c SNe 
come from naked helium stars formed from relatively low-mass progenitors (10 — 2OM0) stripped 
by the interaction, while single SNe Ib/c arise mainly from more massive stars > 24Mq which 
become Wolf-Rayet stars prior to explosion (see Fig. 2a). In a model containing only single stars 
with Z — 0.02 the SN Ib/c rate is reduced to 2 x lO-^yr"^ which would be consistent with ob- 
servations. These rates are sensitive to the assumed mass-loss rate for red supergiants (see Hurley 
et al. 2000), as well as to our assumption that black hole formation ignites a SN explosion, both 
of which are very uncertain. If mass-loss rates are lower or if black holes form without a SN, the 
SNIb/c rate from single stars would be lower. 
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3.2 Gamma-ray bursts 

GRBs are thought to be produced when the evolved core of a massive star collapses to a black 
hole and the remaining star has too much angular momentum to fall in directly. Using the rapid 
binary code we are able to generate a series of large single and binary populations and evaluate 
the formation rate of interesting core collapse species. Rotating naked helium stars, presumed to 
have lost their envelopes to winds or companions, are evolved from the ZAMS ignition upto the 
formation of the final CO core. 

The cumulative rate (above a threshold core mass, Mco) for stars of Z = Zq that explode as 
type Ib/c SNe is shown in Fig. 2a. A fraction, /bh, of the above stars will form a central black 
hole after exploding (see Fig. 2b). By including angular momentum transport by non-magnetic 
processes, such as mass loss, mass transfer and mass accretion, we are able to roughly estimate the 
fraction /q of SN Ib/c progenitors that end their lives with sufficient angular momentum to form a 
centrifugally supported disc (see Fig. 2c). 

High angular momentum is a common requirement in all current GRB models that involve 
massive stars, including those that use a pulsar power source. Yet the actual angular momentum 
in a presupemova star is unknown. Without magnetic fields, and with our approximate treatment 
of angular momentum transport, we find that a naked helium star cannot retain enough angular 
momentum to form a centrifugally supported disc around the collapsed object, as required by the 
coUapsar model of GRBs (fa ^ for single stars; see Fig. 2c). The main reason is that during 
the early phase of WR evolution, the progenitor star spins down significantly as it loses a large 
fraction of its mass quite rapidly (Fig. 3). This is even more so if magnetic fields couple the core 
effectively to the envelope - the magnetic interaction of the rapidly rotating helium core with its 
stationary envelope during the red supergiant phase will halt the core, making it unfit as a coUapsar 
progenitor (Spruit & Phinney 1998; Spruit 2002; Heger & Woosley 2002; but see Livio & Pringle 
1998). Either the description of the stellar model is inaccurate or some route other than single star 
evolution must be involved in making GRBs. The natural alternative is a close binary. 

If either star fills its Roche lobe, then gas flows from the outer layers of the star through the 
inner Lagrangian point that connects the two Roche lobes. Some or all of this gas may be captured 
by the companion star so that mass transfer occurs and, as a result, the subsequent evolution of 
both stars takes a different course from that of isolated stars with important consequences for their 
orbit and spin (Fig. 4). If the hydrogen envelope is removed too early (or the two stars coalescence) 
mass loss from the WR star also reduces the angular momentum (Fig. 5). The effect that the com- 
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panion mass has on the outcome of binary evolution is illustrated by comparing the evolution of 
the primary star in Figs. 4 and 5. 

Finally, the expected birth rates of stars in the Galaxy that have both a massive enough core to 
form a black hole and an adequate rotation rate at the time of collapse to allow the formation of 
a disc are shown in Fig. 2d. These core collapse species may be the progenitors of the common 
(long-soft) GRBs. 

3.3 Sensitivity to model parameters 

To give an idea of how sensitive the final state of the system is to changes in the physical parameters 
that govern the evolution, we reconsider the example illustrated in Fig. 2. If the common-envelope 
efficiency of the component stars is taken to be acE = 3, rather than acE = 1, less energy is 
required to drive off the common envelope. The result is similar to a reduction in the envelope 
binding energy and so it is less likely that the process ends in coalescence of the cores. This 
increment in the common-envelope efficiency parameter slightly increases the GRB production 
rate (see Fig. 6). Wider systems are not brought within an interaction distance after the common- 
envelope phase. Thus systems that do interact when ckce = 1 do not do so when ckce = 3. 
However, the net increase is mainly because closer systems that would coalesce when ccce = 1 
now remain detached. 

At first sight this may seem an unphysical model given the definition but, as discussed by 
Iben & Livio (1993), an increase in acE can be possible if additional energy sources other than 
the orbital energy are involved. Processes with potential to supply such energy include enhanced 
nuclear burning in shell burning zones of giants, nuclear burning on the surface of a degenerate 
secondary, dynamo generation of magnetic fields and recombination of the hydrogen and helium 
ionization zones in giants. Possibly the common envelope absorbs ordinary nuclear energy in the 
process of swelling up, but this should occur on a thermal timescale. Unfortunately, the theoretical 
determination of reliable values for acE has proven difficult owing to a lack of understanding of 
the processes involved and our inability to model them. 

An increase (or decrease) of the mass-loss rate has profound consequences for the WR popu- 
lation. At a given metallicity, the minimum initial mass for the formation of a WR star (lower for 
higher M), the duration of the WR stage (greater for higher M), the times spent in the different 
WR subtypes and the surface composition during these phases are all very sensitive to the mass- 
loss rates (Maeder 1991). There is no doubt that major changes to the outcome result if Mwr is 
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altered by even a small fraction Jm- Inspection of Fig. 7 reveals that as the mass-loss is reduced, 
the threshold for the removal of the hydrogen envelope by stellar winds is raised and the loss of 
angular momentum is inhibited. This increases the mass of the CO core and favours black hole 
formation. Conversely, increasing the mass-loss rate promotes the loss of angular momentum and 
inhibits black hole formation (see Fig. 7). 

It should be noted that the choice of the initial mass above which stars become black holes 
rather than neutron stars is not well constrained. The low-mass X-ray binary (LMXB) catalogue 
compiled by van Paradijs (1995) lists about 80 bright persistent sources in the Galaxy but makes 
no distinction as to whether the compact star is a neutron star or black hole. While this number is 
in fair agreement (HOO) with the birth rates derived with the prescription for the remnant's mass 
given in equation (17), the differences introduced by this choice in the overall evolution are small 
when compared with those resulting from varying the (very uncertain) initial distributions. The 
birth rate of persistent neutron star LMXBs derived by means of such a prescription is similar 
to the rate found by Portegies Zwart et al. (1996) in their standard model that includes velocity 
kicks. Measurements of black hole masses in binaries, although still highly uncertain, yield values 
between 3 and 2OM0 (Orosz et al. 2001; McClintock et al. 2001; Froning & Robinson 2001; 
Wagner et al. 2001; Table 1 of Fryer & Kalogera 2001 and references therein). For progenitors of 
40 ^ Mzams/-^© ^ 100, equation (17) gives final remnant masses in the range 1.8 — S.OM© (Fig. 
8). If the amount of mass that a black hole can accrete is limited by the Eddington limit (Cameron 
& Mock 1967) then a black hole is unlikely to accrete enough material to take its mass above 
5Mq and this model is inconsistent with the observations. There is, however, some uncertainty as 
to whether the Eddington limit should actually be applied because the energy generated in excess 
of the limit might be removed from the system in a strong wind or asymmetrically through a 
disc. Super-Eddington accretion rates may be crucial in allowing black holes in binaries to gain a 
significant amount of mass. 

Imposing the Eddington limit significantly reduces our ability to grow black hole masses to 
the observed values and, for this reason, in addition to our somewhat arbitrary prescription (see 
equation 17) for neutron star and black hole masses, we also use the formalism developed by 
Belczynski, Kalogera & Bulik (2002b; hereafter B02). Based on the Hurley et al. (2000) formulae, 
they calculate the final CO core of a star and then use stellar models of Woosley (1986) to obtain 
a final FeNi core mass as a function of the CO core mass (see their equation [1]). To estimate the 
mass of the remnant formed, B02 follow the results of hydrodynamical core-collapse calculations 
(Fryer et al. 1999; Fryer & Kalogera 2001). Fig. 9 compares the black hole birth rates derived 
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using the B02 formalism with those based on equation (17). The birth rate of black holes shows 
an increase for stars with low CO core masses. Comparing the birth rates of GRBs in both models 
reveals a variation by a factor of 5. The relevance of the relative distribution of neutron stars and 
black holes on the birth rates of GRBs is further complicated by the fact that a rapidly rotating 
neutron stars with an ultrahigh magnetic field could, in principle, also serve as a trigger (Usov 
1992; Thompson 1994). 

These examples serve to demonstrate the sensitivity of binary evolution to the choice of model 
parameters. 

3.4 Stellar companions 

One of the most important questions relating to WR stars is whether they are all members of binary 
systems or rather do some truly single objects exist. Possibly the majority of stars are members 
of binary or multiple systems but most are sufficiently far from their companions that their evo- 
lution proceeds essentially as if they were single stars. However there is an important minority 
of stars which are close enough that their evolution is dramatically changed by the presence of a 
companion. The most dramatic effects occur when one of the stars in the system is a very compact 
object. Here we investigate the expected binary systems with black hole primaries which may be 
the sources of GRBs earlier in their evolution. 

Table 2 shows the various types of companions to the GRB progenitor star at the time of col- 
lapse. For both secondary mass distributions the vast majority of companions are, without surprise, 
MS stars and thus likely donors of soft X-ray binaries (Lee et al. 2002). The typical orbital pe- 
riods for these binaries at the time of explosion are in the 0.01 - 1 days range (Fig. 10). Most 
companions, being MS stars, have small mass-loss rates of lO^^Meyr^^ or less and are thus un- 
likely to significantly enhance the density around the GRB progenitor. There is a small, interesting 
sub-class of binaries in which the companion star is a collapsed object (about 1% for the KTG-q 
distribution). In this case, the GRB, which is likely caused by the relativistic jet, expelled along 
the rotation axis of the collapsing stellar core, propagates within a SN remnant of about 1 Myr in 
age. 

4 THE ROLE OF METALLICITY 

The assumption that all stars within the population are bom with the same composition is some- 
what naive: nucleosynthesis in successive generations of stars enriches the gas from which they 
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form as the Galaxy evolves. Metallicity Z influences the stellar evolution of massive stars mainly 
through bound-free and line opacities in the outer layers of massive stars owing to their influence 
on stellar wind mass loss rates. The numbers of WR stars in galaxies of different metallicities are 
an important test of the stellar models at various Z and of the values of the final stellar masses. The 
theoretical predictions of the ratios of WR/O-stars, WCAVR and WCAVN are in good agreement 
with observations of the Local Group (Maeder 1991), which support the following dependence of 
the mass-loss rates on metallicity: (MJMq) = {Z/Zq)°-^ (de Jager et al. 1988). 

From our large binary population models with a variety of initial masses and metallicities we 
find that both the fraction of stars that form a central black hole and those that end their lives 
with sufficient angular momentum to form a disc significantly increase with decreasing metallicity 
(Figs. 11 - 13). This is because low metallicity keeps the radius of the star smaller and reduces 
mass loss. Both properties inhibit the loss of angular momentum (MacFadyen & Woosley 1999; 
Heger & Woosley 2002) so low-Z stars are likely to be rotating more rapidly. Also, the lower 
the metallicity, the higher the stellar mass for WR star formation. The latter increases the mass 
of the heaviest CO core and favours black hole formation which in turn could power a are more 
luminous burst because more energy can be extracted from the black hole via e.g. the Blandford- 
Znajek mechanism (Blandford & Znajek 1977; Lee et al. 2000; Wheeler, Meier & Wilson 2002). 

The likely metallicity dependence of both black-hole formation and rotation suggests that the 
production of GRBs is likely to affected by the physical conditions in the local ISM - at high Z gas 
opacities are larger in the outer stellar layers and so more momentum is transferred by radiation 
pressure, mass-loss is more intense and GRB formation is disfavoured. 

What are the effects of a dependence of GRB luminosity on the metallicity of their progeni- 
tors? The most significant is a potential offset between the true star-formation rate and that traced 
by GRBs (Ramirez-Ruiz, Lazzati & Blain 2002). If GRBs in low-metallicity environments and in 
low-mass galaxies are more luminous then they are likely to be overrepresented in GRB samples. 
Low-mass galaxies are likely to have statistically lower metallicities and thus contain more lumi- 
nous GRBs than high-mass galaxies. Because galaxy mass is expected to build up monotonically 
through mergers, it is possible that the highest- 2; GRBs could be systematically more luminous 
due to the lower mass of their hosts. This effect is likely to be more significant than, but in the 
same direction as, the global increase in metallicity with cosmic time. 

The most luminous GRBs of all could be associated with metal-free Population-Ill stars but 
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their very high redshifts would make examples difficult to find even in the Swift'^ catalogue of 
hundreds of bursts. 

Star-formation activity is likely to be enhanced in merging galaxies. In major mergers of gas- 
rich spiral galaxies this enhancement takes place primarily in the inner kpc, as bar instabilities 
drive gas into the core (Mihos & Hemquist 1994). Metallicity gradients in the gas are likely to be 
smoothed out both by mixing prior to star formation and by SN enrichment during the burst of 
activity. GRB luminosities could be suppressed in such well-mixed galaxies making GRBs more 
difficult to detect in these most luminous objects. 

Shocks in tidal tails associated with merging galaxies are also likely to precipitate the forma- 
tion of high-mass stars, yet as such tails are likely to consist of relatively low-metallicity gas, it 
is perhaps these less intense sites of star-formation at large distances from galactic radii that are 
more likely to yield detectable GRBs. This might have the unfortunate consequence of making 
GRBs more difficult to use as clean markers of high- 2; star-formation activity. More optimistically, 
the astrophysics of star formation in high-redshift galaxies could be studied using the intrinsic 
properties of a well-selected population of GRBs with deep, resolved host galaxy images. If there 
is a bias towards the discovery of GRBs in low-metallicity regions, then the GRB host galaxy 
luminosity function will be biased to low luminosities. 

5 COMPARISON WITH OBSERVATIONS 

The high Lorentz factors and energies seen in GRBs are consistent with the catastrophic formation 
of a stellar black hole with mass of a few Mq, with about 1 % of the rest mass energy going to 
a relativistic outflow. This could be the extreme example of the asymmetric explosion produced 
by supernovae (Khokhlov et al. 1999) in which, instead of halting at the neutron star stage, the 
collapse continues to the black-hole stage producing an even faster jet in the process (MacFadyen 

6 Woosley 1999). GRBs arising from a small fraction of stars that undergo this type of catastrophic 
energy release are likely to produce coUimated outflows. 

Even if the outflow is not highly coUimated some beaming is expected because energy would 
channeled preferentially along the rotation axis. Also, one would expect baryon contamination to 
be lowest near the axis because angular momentum flings material away from the axis and mate- 
rial with low-angular momentum falls into the black hole. The dynamics, however, are complex. 



^ http://swift.gsfc.nasa.gov 
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While numerical simulations of collapse scenarios can address the fate of the bulk of the mat- 
ter (MacFadyen & Woosley 1999; MacFadyen, Woosley & Heger 2001; Aloy et al. 2000, 2002; 
Zhang et al. 2003), higher resolution simulations of the outer layers of the stellar mantle seem 
to be required because even a very small number of baryons polluting the outflow could severely 
limit the attainable Lorentz factor. The entrained baryonic mass would need to be below IO^^Mq 
to allow these high relativistic expansion speeds. 

Although jets in GRBs were first suggested for GRB 970508 (Waxman, Kulkarni & Frail 
1998), they were widely invoked for GRB 990123 to explain its spectacular energy release. Sub- 
sequent multi wavelength observations of GRBs have been interpreted as evidence for explosions 
with jet-like geometry (Stanek et al. 1999; Harrison et al. 1999; Castro-Tirado et al. 1999). The 
detection of polarization (e.g., Covino et al. 1999; Wijers et al. 1999) gave further credence to the 
jet hypothesis. 

Because conical fireballs are visible to only a fraction of observers, the true GRB rate is 
Rt = < > -Robs where i?obs is the observed GRB rate and < /j^^ > is the mean of 
the beaming fractions. Frail et al. (2001) recently derived < f^^ >^ 500 for a comprehensive 
sample of GRB afterglows with known distances based on observed broadband breaks in their 
light curves^. This implies, within the uncertainties and possible limitations of such a method, that 
only a small fraction of GRBs are visible to a given observer and therefore the true GRB rate is 
several hundred times larger than the observed one: i?obs ~ 10^^i?sNib/c (Schmidt 2001). Fig. 
14 shows our formation rates (as a function of Z) for GRB progenitors, stars that have both a 
massive enough core to form a black hole and adequate rotation at the time of collapse to allow the 
formation of a disc. Clearly these rates are highly sensitive to the assumed black hole formation 
model. Binary models in the KTG-q case are easily capable of supplying a sufficient number 
of progenitors even if the GRB rate is several hundred times larger than the observed rate. A 
discrepancy lies in the birth rates of coUapsars when using an initial KTG-KTG distribution of 
binaries because the theoretical rates in the HOO case are low but not for the B02 case which are 
in good agreement. This can be reconciled by assuming a smaller mass threshold for black hole 
formation (see equation 16). 

3 This estimate assumes that the breaks observed in many GRB afterglow Ughtciu-ves are due to a geometrical beam effect and not to either 
a transition to non-relativistic expansion or an environmental effect such as a sharp density gradient. An important feature produce by this jet 
coUimation is the achromaticity of the afterglow break, which clearly distinguishes it from the steepening that may be produced by the passage of 
a spectral break through the observing band. 



6 CONCLUSIONS 
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The primary purpose of this paper was to use the rapid binary-evolution algorithm to evaluate 
the formation rate of interesting individual species of WR stars that may be the progenitors of 
the common long-soft GRBs. Rotating naked helium stars, presumed to have lost their envelopes 
in winds or to companions, are followed from the ZAMS up to the formation of the CO core. 
Recognizing that the two essential ingredients for the coUapsar model are a sufficiently massive 
core to form a black hole and enough rotation to form a disc, we have studied their effects on the 
results of binary populations synthesis and made a global comparison with GRB observations. 

The framework we have used to determine a star's angular momentum is a simple one. While 
many of the processes involved are not understood in detail, we do have a qualitative picture of 
how binaries evolve and we hope to construct a model that correctly follows them through the 
various phases of evolution. This method has led us to conclude that the resulting spin rates for 
single stars become too low to form a centrifugally supported disc to drive a GRB engine. Heger 
& Woosley (2002), including more realistic estimates of angular momentum transport, argue that 
some single stars may still be able to form a disc. However inclusion of magnetic torques in their 
calculations also results in too little angular momentum for coUapsars. Either the stellar model 
description is inaccurate or some other evolutionary path must be involved in making GRBs. The 
obvious option is a close binary or merger (see also Fryer et al. 1999), where tidal interaction 
transfers orbital angular momentum to the stellar rotation. An important ingredient is the fact that 
the rapidly rotating new black holes that power the GRBs are typically accompanied by MS stars 
and so are likely donors of soft X-ray binaries (only those MS stars which remained gravitationally 
bound after the explosion). 

For most of the individual binary populations we find that the expected formation rates of 
coUapsars in the Galaxy are easily capable of supplying a sufficient number of progenitors, even 
if the true GRB rate is several hundred times larger than the observed rate because of beaming. 
An interesting conclusion of this work is that binary stars at low metallicity are important for 
the formation of a rapidly rotating, massive helium core at collapse. This effect could increase 
the GRB formation rates by a factor of 5-7 at Z = Zq/ 200. We finally note that the simple 
binary-evolution algorithm provides an adequate description of the observations although more 
stringent constraints to many of the evolution variables are needed in order to draw more accurate 
conclusions. 
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Table 1. Model rates of type Ib/c and type II SNe (per Galaxy per year) for single stars and binaries with the secondary chosen either from the KTG 
or q IMF. SN rates include binaries which may have coalesced or broken up. The associated uncertainties in the observed offsets represent the la 
confidence region as calculated from the Poisson errors. S = 7.608 Galaxy-^ yr-^ is the rate of star formation. 
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Table 2. Stellar type of the companion star when the primary explodes as a possible GRB (percentages) for Z = Zq in the binary phase. 
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Figure 1. Selected evolutionary tracks in the HR diagram of massive stars (M ^ IOMq) with initial Z = 0.02 composition. The different 
line styles correspond to different evolutionary phases: main sequence (MS); Hertzsprung gap (HG); first giant branch (GB); core helium burning 
(CHeB); asymptotic giant branch (AGE); and naked helium star (He), (a) Single stars evolved with mass loss by stellar winds, (b) Binary evolution 
with a secondary IOMq star initially orbiting the primary with a 100 days period. 
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Figure 2. Formation rates of interesting core collapse species, (a) Cumulative rate (above a threshold core mass Mco) for Z = Z© stars that 
explode as type Ib/c SNe. The shaded region represents the la uncertainties on the birth rates of type Ib/c SNe as derived by Cappellaro (2001). 
The birth rates reported above assume Ho = 65 km s~^ Mpc~^, and that the Galaxy has an average SN rate for its morphological type (here 
assumed Sb-Sbc) and luminosity (2.3 X lO^'^Lg q). The thin solid line assumes a 50% binary fraction (KTG-q case), (b) Fraction of SN Ib/c 
& Mco) that form a central black hole after exploding, (c) Fraction of SN Ib/c (^ Mco) th^t end their lives with sufficient angular momentum 
to form a centrifugally supported disc, (d) Cumulative rate (^ Mco) of stars - thought to be GRB progenitors - that have both a massive enough 
core to form a black hole and adequate rotation rate at the time of collapse to allow the formation of a disc. 
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Figure 3. Single stellar evolution after core hydrogen burning for Z = Zq and initial masses 40, 60, 80 Mq . Both the mass and angular momentum 
are shown for the post-MS phase as a function of time. 
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Figure 4. The evolution of the mass and angular momentum of various primary stars after core hydrogen burning for Z 
20, 40, 60 Mq with a binary companion with initial P=100 days and M2 = IOMq. 
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Figure 5. Similar to Fig. 4 but with M2 = 3Mq . 
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Figure 8. The probability distribution of final remnant masses for an initial KTG-q distribution of binaries and for two different prescriptions for 
NS (thick Unes) and BH (thin Unes) masses. See text for details. 
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Figure 9. Similar to Fig. 2 but using the B02 prescription for NS and BH masses ratlier than the HOO formulae. An initial KTG-q distribution of 
binaries is adopted. 
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Figure 10. The probability distribution of binary period at the time when one of the stars explodes. The initial explosion then produces a BH and 
sufficient angular momentum to form a disc. 
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Figure 11. Similar to Fig. 2 for various Z using an initial KTG-q distribution of binaries. 
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Figure 12. Similar to Fig. 1 1 but with an initial KTG-KTG distribution of binaries. 
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Figure 13. Similar to Fig. 1 1 but using the B02 prescription for NS and BH masses rather than the HOO formulae. 
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Figure 14. GRB formation rates as a function of metallicity for two different prescriptions for NS and BH masses. A GRB progenitor is thought to 
be a star that has both a massive enough core to form a BH and adequate rotation rate at the time of collapse to allow the formation of a centrifugaUy 
supported disc (see Section 3). The shaded region marks the la width of the birth rate of GRBs as derived by Frail et al. (2001). The filled circles 
are KTG-q binary populations, the empty circles are binaries with an initial KTG-KTG distribution. 



